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Abstract. In this work we consider a DBI Galileon (DBIG) inflationary model and constrain 
its parameter space with the Planck 2015 and BICEP2/Keck array and Planck (BKP) joint 
analysis data by means of a potential independent analysis. We focus our attention on 
inflationary solutions characterized by a constant or varying sound speed as well as warp 
factor. We impose bounds on stringy aspects of the model, such as the warp factor (/) 
and the induced gravity parameter (m). We study the parameter space of the model and 
find that the tensor-to-scalar ratio can be as low as r ~ 6 x 10“^ and inflation happens 
to be at GUT scale. In addition, we obtain the tilt of the tensor power spectrum and test 
the standard inflationary consistency relation (r = —8nt) against the latest bounds from the 
combined results of BKP+Laser Interferometer Gravitational-Waves Observatory (LIGO), 
and find that DBIG inflation predicts a red spectral index for the tensor power spectrum. 
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1 Introduction 

In the context of the recent results from Planck [1-3] , the joint analysis of BICEP2/Keck Array 
and Planck (BKP) [4], the inflationary paradigm, in particular single-field inflation, seems 
to be the one chosen by nature to generate the observed adiabatic, nearly scale-invariant, 
Gaussian spectrum of curvature perturbations. In particular, the Planck satellite, with its 
exceptional quality data, constrains the power spectrum tilt of the curvature perturbation 
with Us = 0.968 ± 0.006, ruling out scale invariance at more than 5 ct [2]. Moreover, re¬ 
sults indicate that we live in a spatially flat universe, Dfc = 0.000 ± 0.0025, and that the 
perturbation spectrum imprinted in the GMB is Gaussian to a high degree. This imposes 
severe constraints on the bispectrum amplitudes: = 0.8 ± 5.0, = —4 ± 43 and 

fN^° = —26 ± 21 at 68% confidence level. Furthermore, Planck data [2] suggests a small 
running dns/dlnk = —0.003±0.007 , which is consistent with the prediction from single-field 
models of inflation [5]. Although there is no significant detection of primordial tensor modes, 
the current analysis points to an upper bound for tensor-to-scalar ratio (r) of the order of 
O (l0“^). The proposed post-Planck satellites GMBPol, GOrE, Prism, LiteBIRD and many 
other ground based experiments such as Keck/BIGEP3 [6, 7] are expected to reach enough 
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sensitivity to detect B-modes and further constrain r ^ O (l0“^). The latest results sug¬ 
gest no evidence for a blue tilt of the gravitational wave power spectra [2, 4], In addition, 
this result is supported by the Joint analysis of BKP+Laser Interferometer Gravitational- 
Waves Observatory (LIGO) [8]. According to the new data, the simplest model with a (fP' 
potential is ruled out, whereas Starobinsky’s Bp inflation and non-minimal coupled models 
[2, 9, 10] remain consistent. In recent works, these two models have been presented as cos¬ 
mological attractor models [11]. Under the new constraints on inflation with a 

non-canonical kinetic term has gained importance, in particular the Dirac-Born-Infeld (DBI) 
inflation is compatible with Planck data when the sound speed Cg > 0.087 [1]. This is in¬ 
teresting in the context of consistent inflationary mechanism aiming UV-completion theories, 
such as string theory and supergravity (SUGRA) [12-15]. 

Our paper provides an observationally consistent inflationary scenario bearing features 
of modified gravity and string theory. We focus on a particular class of string inflation 
constituting a generic extension of the DBI inflation [16-24]. A peculiar feature of the DBI 
action is that it involves a non-linear function of the inflaton’s kinetic term, which leads to 
the generation of a potentially measurable non-Gaussian signal in the perturbation spectrum. 
This is due to a small sound speed which can accommodate a smaller tensor-to-scalar ratio 
while sourcing large non-Gaussianities. However, stringent observational bounds on non- 
Gaussianity restrict the validity domain of DBI inflation [1, 13, 21]. In this paper, we focus 
on a well-motivated extension of DBI inflation: the DBI Galileon (DBIG) inflation [25-27]. 
This model is more generic in the sense that it involves induced gravity, which arises due to 
the motion of the D-brane in warped geometry. 

DBIG inflation, as it has been generically considered in the literature [28-35], is mainly 
restricted to the slow-roll regime. Moreover, these studies are more focussed on the parameter 
space of the single-field and multifield DBIG model with respect to the various types of non- 
Gaussianities. We must emphasize that in Ref. [32] a different motivation is considered, and 
DBIG inflation is studied in the background of SUGRA under the assumption of a Goleman- 
Weinberg type of Galileon potential. In this paper, we propose to study single-field DBIG 
inflation without any particular choice of potential. More precisely, we aim to constrain the 
parameter space of the DBIG model according to the latest Planck 2015 data. We mainly 
focus our attention in two inflationary regimes. Namely, those with and without a constant 
warp factor. We aim to identify crucial differences between these two scenarios with respect to 
the corresponding inflationary predictions. In addition, in each case, we analyze the deviation 
from the standard slow-roll consistency relation r = —8nt due to the effect of induced gravity 
on the D-brane. 

The organization of this paper is as follows. In Sec. 2 we briefly describe the model and 
present the background equations for the DBIG inflation with non-trivial warping [31]. In 
the case of constant sound speed and warp factor, we obtain the exact background solutions. 
In Sec. 3 we compute the amplitude and tilt of the scalar and tensor perturbation spectra for 
the DBIG model, deferring the details of the computation to Appendix B. In Sec. 4 we study 
the parameter space of the DBIG model by comparing its predictions in different limits with 
GMB data. In Sec. 5 we present general background solutions using two different ansatz to 
integrate analytically the equations of motion. A detailed computation of the approximate 
solutions can be found in Appendix A. Finally, we present our conclusions in Sec. 6. 
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2 DBI-Galileon inflationary model 


We begin by reviewing the DBIG inflationary scenario following Ref. [31]. Such a setup 
considers a D3-brane with tension T 3 evolving in a ten dimensional geometry described by 
the metric, 


ds"^ = h {y^) g^udx^dx'" + {y^) Gjj (y^) dy^dy^ = HABdY^dY 


B 


( 2 . 1 ) 


with coordinates Y^ = where // = 0,....3 and I = 1,....,6. The induced metric on 

the D3-brane is given by 

l^u = HABd^Y^%d,Y^%, ( 2 . 2 ) 

where the brane is embedded in higher dimensions by means of the functions Y^-^ with the 
x^ being the space time coordinates on the brane. In brane inflation, the role of the inflaton is 
played by the radial coordinate (p) of the brane that is moving in the extra dimensions. Since 
we are only considering single-field inflation in this paper, we choose the brane embedding as 
'^{b) ~ ^ (x^)). Then the induced metric can be written as 

{g^lu + fd^yduy) , 

where / and (p are the warp factor and the scalar held defined by 

h 


(2.3) 




(2.4) 


The D3-brane here is embedded in 5D geometry with the induced metric Eq.(2.3). This 
introduces an additional contribution in the action known as Galileon term [36]. The total 
action is then given by 


S= d^x 


b] + [7] + V^gJ^br 


(2.5) 


where mp = 2.24 x 10^® GeV is the four dimensional reduced Planck mass, m is a parameter 
associated with the induced gravity^ and 


d^brane — 


f{<P) 




where 


V = det ((5(( -k fd^pdup) ■ 

Assuming the flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric 

ds'^ = —dt^ + a^{t)dx‘^ . 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


and allowing the warp factor / to vary, the gravitational held equations for the action in 
Fq. (2.5) are [31] 


+ = =-[ — -1] + V 


, m 


1 / 1 




/ \cv 


(2.9) 


non trivially depends on the warping h, see [31]. In this paper, rh is treated as a model parameter. 
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— rupH + 


cv 


H 


Cv 


'/iV4' 


/iV4 I cv 


H 3 


1 - 1 ^^ 


a 


( 2 . 10 ) 


where c|, = 1 — is the squared sound speed^, H = H — ^ and . The 

appearance of Eq. (2.10) can be simplihed to 


H - XiH^ + A 2 = 0 

after introducing the functions 

^fiVf - e) 


( 2 . 11 ) 


Ai = 


A 2 = 


m 


rripCv + 

1 - c|, 


din o 


-4 +77 Hr-1 1-4 


2 V c2 

\Cp 


4 


2/ [m?pCv + H) ’ 

which depend on fn and cv- We also introduce the slow-roll parameters 


, ( 2 . 12 ) 

(2.13) 


H 


dine 


^ H ^ dlno 


, Cv = 


dlncx) 
din a 


, m = 


din ev 
din a 


4 = 


dlnf 

din a 


r]f = 


dlnef 


din a 


(2.14) 


to describe the evolution of the background geometry, the sound speed and the warp factor. 
Note also that in the above we take the brane tension T 3 to be a constant, as is usually 
considered. 

In the following we obtain solutions to the background equations for the cases when Ai ^2 
are constants. 

2.1 Constant sound speed and warp factor 

Whenever the sound speed {cv < 1) and the warp factor are constant, i.e., ev = Cf = 0, the 
coefficients Ai ^2 in Eq. (2.11) are constants. Integrating Eq. (2.11) in that case is straightfor¬ 
ward. We obtain 

H = ^ + (2.15) 

Ai 

where k 7 ^ 0 is an arbitrary, dimensionful constant. Writing H = d/a, the solution to 
Eq. (2.15) is 


a2+(t) = 


A 2 


4il^l 

where we introduce 


exp [41 + < 71 ) 7r/2] sec4 {t — — i {1 + (Ti) vjA , (2.16) 


fJi = sign(K) = sign(id) , fT 2 = sign(d). (2-17) 

The explicit time-dependence of the Hubble parameter can be obtained from Eq. (2.16) 


'A 

dd(t) = — ( 4 j (j 2 tanh 


\/AiA2(T2 (t -t) - i{l + fJi) 7r/4 


(2.18) 


^Note that the sound speed ct> here depends not only on the brane dynamics, (as in DBI models [18, 19, 21]) 
but also on the induced gravity [31]. 
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To study inflation we need to set <72 = sign (d) = +1, regardless of ui = sign(ilr). An 
increasing expansion rate is obtained for a\ = +1 (A 2 < XiH^), which corresponds to the 
singular behaviour of the scale factor and the Hubble parameter at t —>• t (purple line) 
displayed in Fig. 1. A decreasing expansion rate corresponds to cJi = —1 (A 2 > XiH^)^ in 
which case both a{t) and H{t) remain finite throughout the entire evolution (blue line). In 
the context of inflation, we focus only on the decreasing expansion rate ai = —1, for which 
we find a non-singular behaviour for the scale factor a{t). 




V ''■2 0-2(t-Tj 


^1 ^2 (f-f) 


Figure 1. Evolution of the scale factor according to Eq. (2.16) (left panel) and the Hubble parameter 
H, according to Eq. (2.18) (right panel). 


In Sec. 4.1 we impose the necessary conditions to obtain an inflationary expansion in 
agreement with current observations. To do so, in the next section we investigate the scalar 
and tensor perturbation spectra, which depend on the slow-roll parameters e and t]. Using 
Eq. (2.14) and Eq. (2.18) we obtain 


e{t) = Aicsch^ \/ XiX2<72{t — t) — z(l -|- (Ti)7r/4 
r]{t) = 2Aicoth^ \/AiA 2 cr 2 (t — t) — i{l -|- cJi)7r/4 
from which we arrive at the relations 


77 = 2(e + Ai) , H^ = X2{Xi + e) 


-1 


(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 


where we emphasize that the slow-roll parameter rj explicitly depends on Ai. During inflation, 
r/ <C 1 implies Ai <C 1. Therefore, several constraints (to be discussed later on) must be 
imposed on the model parameters to have Ai <C 1. 


3 Perturbation spectra 

In this section we portray the generalized approach for calculating the scalar and tensor 
power spectrum as described in the context of generalized G-inflation [37]. We present the 
calculations of spectral index (n^), tensor-to-scalar ratio (r) and tensor tilt (rit) up to the 
slow-roll approximation. 


3.1 Scalar spectrum 


The second order action for scalar perturbations is given by [31, 37], 



dt d?x [a^Gs 



J^s/Gs 


(VC)^ 
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(3.1) 



















where ( is the curvature perturbation, t is the cosmic time, J^s, Gs are arbitrary functions of 
time and Cs = is the sound speed for scalar perturbations. For the DBIG model 

[31], the functions J^s and Qg that determine the second order action for scalar perturbations 
are 


m 


2 r 


J^s{cv, ev, e) = B{t) = nip (e/C(3/C - 2 ) + /C - 1 ) H- 

cv 


(e + eT>)JC f “2 - 2 ] + /C — Cp 


V c 2 


Gs(ct>, ex>) e) = ^(j') — + 3c|,(l — H—g- 


2 r 




-D L 


(e + eD)}C^ + 3c^ ( 1 -^ 


2 \ 1 


(3.2) 


1 mp+c„ 

where H, = — 33 — 5 -. 

rrip+Cjj m 

In addition to the slow-roll parameters previously dehned in Eq. (2.14), we introduce 
the following new parameters 


^ dlnJF, (2) ^ dlnfs 

^ din a ’ * dIn a 

Moreover, using the dehnition of Cg we have 


_ dine* _ 1 ^ ^ 

- ,, — r, \Js 9s) 

din a 2 


9s = 


dine. 


din Qg 
din a 


( 2 ) ^ dlngg 
d In a 


1 


9s 


din a 2 e, 


(/./i"> 


9s9. 


( 2 ) 


(3.3) 


(3.4) 


Although the procedure to obtain the perturbation spectrum is well known, for the sake of 
completeness we include in Appendix B a detailed computation, obtaining [37] 




7s Gsi 


_^ - g2^,-3 r {l^s) ^ 

2 jp3/2 47^2 ’ r(3/2)2 


1 I 9s* fs* 


(3.5) 


where labels the time of sound horizon crossing when ki/g = — 1. If slow-roll parameters 
are sufficiently small one can consider the linear approximation [38, 39], and then the scalar 
spectral index in such case is given by 


n, - 1 ~ 


4e^ 4“ l^fs* 9s* 


(3.6) 


2 -\- 2 e^ fg:^ Qg^ 

3.2 Tensor spectrum 

Similarly to the case of scalar perturbations, the second order action for tensor perturbations 
can be written as [37] 


Sf > = 1 


dt d^x [a^Gt) 


h%- 


Tt/Gt 


(V/i, 




(3.7) 


where Tt and Qt are functions of time and q = {J-tlGtY^"^ is the sound speed for tensor 
perturbations. In the DBIG model, the second order action is determined by the functions 


m 


Tt{cv) = rrip + m^cv , Gt{cv) = m-p H-. 

Cv 

Similarly to Eq. (3.3), we consider now the additional slow-roll parameters 




din Ft 
din a 


(2) ^ din ft 
* din a 


9t = 


d In Gt 

din a 


( 2 ) _ dln^t 
“ dlna 


(3.8) 


(3.9) 
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and using the definition of ct we also have 


_ dinct _ 1 

— ,, — „ Ut 9t) 

dlno 2 


m = 


din et 


t 


d In a 2et 

It is easy to show now that the tensor spectrum is (see Appendix B) 


(2) n 

t - 9tgt 


(^ 1/2 7^2 

* -r3/2 47^2 

•'t* 


^ r(3/2)2 


-1 I 


(3.10) 


(3.11) 


If slow-roll parameters are sufficiently small one can consider the linear approximation and 
obtain the tensor tilt as [38, 39] 


4e* + 3 ft* — 9t* 

—2 + 2e* + ft* — 9t* 


(3.12) 


where the subindex indicates the time of sound horizon crossing, determined by the 
condition kyt = — 1. 

Finally, to further constrain the model parameters with observations we obtain the 
tensor-to-scalar ratio 


r = 


Vt 74 

— = 16 — 
'Pc Is 



(3.13) 


From Eqs. (3.12) and (3.13) we observe that the standard consistency relation of single-field 
inflation, r = —8nt, is in general violated. In Ref. [40] it has been shown that, in the case 
of power law G-inflation, one can have either r > —Srit or r < —8nt depending on the 
model parameters. However, the requirement of subluminal propagation speed of the scalar 
perturbations restricts r < —^nt- 


4 Comparison to observations 

In this section we study in detail the observational predictions of DBIG inflation and examine 
the status of the tensor consistency relation. We study the different limits of DBIG inflation 
beyond the slow-roll approximation and evaluate the effect of higher order corrections in 
slow-roll parameters on the model predictions. 

We explore the parameter space (c® , m, f) of DBIG inflation using the Planck con¬ 
straints on {us, r) and the observed amplitude of the power spectrum Vq ^ ~ 2.2 X 10 ® at the 
pivot scale k* = 0.002 Mpc“^ [2]. In all cases, we find that the predictions of (n^, r) do not 
explicitly depend on the warp factor. Therefore, we first find the range of model parameters 
(cx), m) compatible with the observed values of Ug = 0.968 ± 0.006 and r < 0.1 at the 95% 
GL [2]. After that, we calculate the tensor tilt (nt) for the same parameter space that was 
previously constrained. We expect to find departures from the consistency relation of single- 
field inflation, r = —8nt. Finally, we compare our results with the BKP+LIGO constraints 
on the tensor tilt nt = ~OP^-o '52 Q8% GL [2, 8]. 

4.1 Constant sound speed and warp factor 

Let us examine the parameter space of DBIG inflation with ej) = e/ = 0 in different limits. 
For this we use the solutions derived in Sec. 2.1. We focus only on the decreasing expansion 
rate ai = —1, for which we find a non-singular behaviour for the scale factor a{t). 
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Firstly, the number of e-foldings during inflation can be computed as 





(4.1) 


where t* is the time when cosmological scales exit the horizon and te signals the end of 
inflation, set through the condition e(te) = 1- According to observations, the length of the 
inflationary phase required to solve the flatness and horizon problems is around A^* = 40 to 
= 70. Using Eq. (2.18) and the condition e = 1 to determine te, we integrate Eq. (4.1) to 
obtain 


1 cosh\^/X[X^a2iU - t) - i{l + ai)7r/4] 

N^ = — In -^^ 

Al y/l Al 


(4.2) 


which we can relate to the slow-roll parameters e and 77 = 2 (e -|- Ai) at the time of horizon 
crossing 


_Al_ 

(1 -I- Al) exp[2AiA^*] — 1 


(4.3) 


Using Eqs. (2.9) and (2.21), we hnd the scalar potential V in terms of the model parameters 


Ai + e 


rup + 



1 

7 



(4.4) 


which allows us to hnd the energy scale of inhation after evaluating at the time of horizon 
crossing for cosmological scales. Also, we obtain the mass squared of the inhaton 


where 




(4.5) 


(4.6) 


4.1.1 DBI limit: m —)> 0 

The phenomenology of DBI inhation has been done in recent literature [41, 42] assuming a 
particular form of potential. We emphasize here, however, that in our study we do not assume 
any form of the potential. 

In this limit Ai —?■ 0 (see Eq. (2.12)), and we obtain the corresponding background 
solution from the one obtained in Sec. 2.1 as the zeroth order in a series expansion around 
Al = 0. Operating similarly for the number of e-foldings in Eq. (4.2) we easily obtain 


1 - 4 

2fm?pCv 




A 2 


1 


1 + 2N 




where N is the number of e-foldings before the end of inhation. Fixing the number of e-foldings 
and the amplitude of the perturbation spectrum we constrain the warp factor /. Since we 
treat cp as a model parameter, we obtain its range from the prediction for non-Gaussianity 
( 73 — models [16, 17]. Although more accurate expressions exist in the 

literature [28-31, 33, 34, 43], for our purposes it suffices to consider this simple estimate. This 
is appropriate since in the absence of a clear detection of non-Gaussianity [1], the use of more 
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elaborate or complicated expressions is, in principle, uncalled for. Therefore, in this paper we 
will not be concerned with non-Gaussian computations and will use the above expression to 
constrain the sound speed c®. The analysis of the Planck data on r < 0.1 and = —4±43 
allows to set a conservative bound for this 0.087 < cx) < 0.6 [1, 2]. Note that larger values 
of cx>, albeit allowed by the bound from non-Gaussianity, are disfavoured as they result in a 
tensor-to-scalar ratio in excess of the current bound r < 0.1. Fig. 2 represents the viability 
of the DBI model. Because of the stringent bound on the DBI inflation is not capable 
to induce r < 0.01 which is consistent with previous studies [21, 44]. The range of model 
parameters obtained for 0.087 < cx) < 0.6 can be found in Table 1. In Fig. 2 we depict our 
results in the DBI limit. 


0.15 


0.10 


0.05 


0 . 


8%50 0.955 0.960 0.965 0.970 0.975 0.980 

ris 



Figure 2. In the left panel we depict tensor-to-scalar ratio vs. spectral index where in the plot Ah 
varies from 50 to 60 (from left to right) and cx> varies from 0.087 to 0.6 (from bottom to top). In the 
right panel we plot the ratio r/rit vs. sound speed c-p for — 60 . 


4.1.2 Galileon limit: fh^mp 


Although studying this limit is not generic with respect to the structure of DBIG, 
nevertheless be useful to understand the role of induced gravity. Since cx) < 1, 
gives 



A2 = 


1 - 

2cvfm^ ' 


this would 
Eq. (2.11) 

(4.8) 


The slow-roll parameters in this case which are given below 


e 


3(1-cl) 



- 2c 


2 

V 



-\- 2e. 


(4.9) 


Unlike in the DBI limit (cf. Eq. (4.7)), in the Galileon limit, the slow-roll parameters 
explicitly depend on the sound speed. It is obvious from Eq. (4.9) that cx) 1 would 
actually spoil the smallness of rj. Therefore, in this case we need to keep the sound speed 
in the narrow range 0.995 < cp < 1 for the results to agree with the current Planck data. 
Any value of cp < 0.995 would essentially spoil the prediction of the spectral index and 
its value would be signihcantly out of the current bounds Us = 0.968 ± 0.006. Therefore, 
observationally viable inflation due to the induced gravity term sets cx) < 1, thus resulting in 
small non-Gaussianities. This allows to discriminate between the current case and the DBI 
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limit previously studied. Also, the consistency of the predictions with data becomes better 
as the number of e-foldings reduces. In particular, for ~ 50 our results are perfectly 
consistent with current data whereas for > 60 the model is ruled out. Our results in this 
case are depicted in Fig. 3. The derived model parameters for 0.995 < Cd < 1 can be found 
in Table 1. 



(192 0.93 0.94 0.95 0.96 0.97 0.98 0.99 1.00 


ly 



CD 


Figure 3. Plots of spectral index rig vs. tensor-to-scalar ratio r (left) and the ratio r/rit vs. sound 
speed CD (right) in the Galileon limit. In the left panel, we take A* varying from 50 to 60 (from 
bottom to top). For the right panel we considered A* = 60. 


4.1.3 DBI-Galileon case 

In this section we consider Einstein and Galileon gravity are on an equal footing. In this case 


Ai = 


m 


m?pCv + 


-v 


A 2 = 


1 - 


2 / {m?pcv + m^) 


(4.10) 


The corresponding slow-roll parameters are (expressing in the units of mp = 1) 


3 (1 — Cp) fh? 

-w- , V 

-^-1 
^T) J 

[2c|, — [cp — 3 ) 771^] e _ 2cp {cp + m?) 

Similarly to the Galileon limit studied in Sec. 4.1.2, the sound speed needs to be tuned 
to cp ~ 0.98 — 0.99 to keep the slow-roll parameter r/ small enough to have Ug = 0.968 ±0.006. 
We find that cp < 0.98 would essentially spoil the prediction of scalar tilt. We also note here 
that Cp = 1 we obtain exact scale invariance, i.e. rig = 1. Since the slow-roll parameter e in 
Eq. (4.11) depends on the parameter m, the tensor-to-scalar ratio varies for different values 
of the induced gravity parameter fh. This allows us to identify the range of the parameters 
consistent with current data. In Fig. 4 we study the parameter space {cp , fh) using the 
bounds on (n^ , r). The plot shows that, in the limit m —)• 0, the model reduces to DBI case. 
Moreover, unless fh < mp, the effect of the induced gravity forces us to constrain the sound 
speed to cx) ~ 1 in order to maintain the agreement with observations. 

To constrain the model parameters {cp , fh) with the bounds of (n^ , r) it is also necessary 
to check if non-Gaussianities are large. Since the full study of non-Gaussianity is beyond the 
scope of this paper, we use the results in Ref. [31], where the authors study non-Gaussianity 
in the multifield DBIG inflation model. We adopt their expression for in the single-field 


(i - 0 

{cp ± m^) 


±2e. (4.11) 
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Figure 4. Contour plots in the plane {rh,CT>) (with m in units of mp). Blue and orange regions 
represent the space where Ug = 0.968 ± 0.006 and 0.01 < r < 0.1, respectively. 


limit, i.e. taking the adiabatic and isocurvature mode transfer function T„s —t 0. We thus 
constrain our parameter space using the approximate expression [31] 

eg _ 5 21 - 404a + 2233a2 - 3066a3 _ fH'^rh^ 

324c|, (1 — 5a)^(l — 9a) ’ c|, 

Setting W = 60, in Fig. 5 we plot the model predictions in the plane (n^, r) (left panel) 
for different values of Cx> and for different ranges of m, as indicated. In the plotted curves, 
the tensor-to-scalar ratio decreases as we increase rh. Therefore, our results show that an 
increase of the induced gravity lowers the tensor-to-scalar ratio. In the right panel we plot 
the ratio r/rit as a function of rh. In the range of values of c® consistent with the observed 
value of the spectral index we find a slight deviation from the standard consistency relation. 
Nevertheless, such a deviation does not seem to be sufficiently significant to be detected with 
confidence. 

In Fig. 6 we plot the mass squared of the inflaton, as obtained from Eq. (4.5) evaluated 
at the time of horizon crossing for cosmological scales (left panel), and calculated from 
Eq. (4.12) (right panel). From the left plot, we find that the inffaton is tachyonic, whereas 
for smaller values of m, we recover a potential with positive curvature, in agreement with the 
DBI case. In this sense, it may be worth mentioning that the authors in Ref. [45] have studied 
the possibility that the Born-Infeld tachyon be equivalent to a scalar field in an effective field 
theory in different warped geometries. Moreover, in Ref. [42] the observational constraints on 
tachyon and DBI inflation were studied, and the authors showed that tachyon inflation fits 
better with cosmological data than DBI. It is also important to notice that < 0 in all cases, 
which is statistically preferred by data after the Planck and BKP joint analysis [2, 4], Also, 
the joint analysis of BKP+LIGO indicates a red tensor tilt nt = ~0.76]'^q 52 at the 68% CL 
[8]. In Table 1 we report the values of the ratio r/m, which only results in a slight deviation 
from the standard consistency relation in most of the cases. We recall that future cosmology 
probes will be able to discriminate inflationary models by direct detection of primordial B- 
modes [7]. Finally, from the right panel of Fig. 6 we find that the non-Gaussianity parameter 
is consistent with the stringent bounds imposed by Planck data [1]. 
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Figure 5. Plots of spectral index Ug vs. tensor-to-scalar ratio r (left panel) and the ratio r/rit vs. 
m (with 771 in units of mp) (right panel) in the DBIG model. In the left panel we take c-d = 0.98 
and 0.3 < m/mp < 0.72 (red), c-d = 0.985 and 0.5 < m/mp < 1.25 (black), cp = 0.99 and 
0.5 < fh/mp < 1.25 (blue). In the plotted curves fh increases as r decreases. In the right panel, the 
plotted curves correspond to c-p = 0.98 (red), cp = 0.985 (black) and cp = 0.99 (blue). 




Figure 6. Plots of the mass squared of the inflaton field (left panel) and the non-Gaussian parameter 
(right panel) as a function of m (with fh in units of mp). In this plot 0.22 < a < 0.32 for 
0.5 < fh < 1.25. We take cp = 0.985 to build the plots, hence the depicted behaviour corresponds to 
the black line in Fig. 5. 


4.2 Varying both sound speed and warp factor 

The cases considered in Sec. 4.1 (constant sound speed and constant warp factor) are consis¬ 
tent with observational data. However, it is interesting to understand the cases with varying 
Cp and /. The questions we can pose in these cases are, can we get a parameter space with 
r ~ O (l0“^)? How do the warped geometries and the scale of inflation change when (cp, /) 
change with time? What is the nature of inflaton field is such cases? In this section, we 
obtain exact background solutions in two cases: a slowly varying sound speed at fixed warp 
factor and a slowly varying warp factor at fixed sound speed. 

4.2.1 Varying sound speed (ep / 0,r7p = 0) and constant warp factor (e/ = 0) 

We assume a slow variation of the sound speed, i.e. ep <C 1. Using the definition of slow-roll 
parameters from Eq. (2.14), we can approximate cp in terms of V = in a as 

Cp = Crfexp(epiV) ~ Crf(l-hepiV), (4.13) 
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where q is a constant whose magnitude is set some four e-foldings after the largest cosmo¬ 
logical scales exit the horizon. 

To integrate the background Eq. (2.11) it is now convenient to rewrite it as 

H' -XiH+^ = ^, (4.14) 

ri 

where Ai ^2 are computed using the approximation in Eq. (4.13) and the prime stands for 
' = Integrating Eq. (4.14) we obtain the solution H = H{N). To hx the integration 
constant in the solution it suffices to impose that e = — ^ = 1 at the end of inffation. 
We choose not to include here the solution H = H{N) as it is a complicated expression 
involving imaginary error functions [46]. To constrain the model parameters we proceed 
as in Sec. 4.1. Since in this case {us , t) do not depend on warp factor /, we may hnd 
the range for (cd, m, e®) using the current bounds on (n^ , r). Since we assume a slowly 
varying sound speed, its constraint in this case is not signihcantly different from the one 
obtained in Sec. 4.1.3. Consequently, we must tune Cd — 0.98 so that the spectral index 
agrees with observations. We also hnd that consistency with observations demands ex) < 0. 
This resembles the result of Ref. [38], where it was shown that DBI inhation with a decreasing 
sound speed results in an expanding universe, in contrast to the case of increasing sound speed. 
The observables in this case (n* , r) are not very different from those obtained for a constant 
sound speed and warp factor in Sec. 4.1.3. In fact, after an extensive numerical study we hnd 
it difficult to obtain r O (10-3) in this case. Therefore, from our analysis we conclude that 
DBIG inhation with a varying sound speed and constant warp factor does not bring any new 
features. 


4.2.2 Varying warp factor (ej ^ 0,r]f = 0) and constant sound speed {ex> = 0) 

In general, the warp factor can depend on helds not stabilised during inhation. Therefore, it 
is feasible to expect a time-dependent warp factor while cosmological scales are exiting the 
horizon. For example, in Ref. [47], various solutions for warped geometries were considered 
in the context of DBI inhation. In the following, we consider a slowly varying warp factor in 
the DBI-Galileon inhation model and constrain its variation using current data. Therefore, 
taking e/ <C 1 we approximate the warp factor as follows 

/ =/oexp(e/iV) ~/o (1-ke/V) , (4.15) 


where /q is the initial value warp factor and ej is constant and treated as free parameter. 
Similarly to the previous case, we set the magnitude of /o four e-foldings after the largest 
cosmological scales exit the horizon. 

It is important to remark that, in contrast to the previous case, where Ai ,2 = Ai, 2 (-V) 
and no simple analytical solution can be found for Eq. (4.14), using €x> = 0 and e/ = const, 
gives Ai = const, and only A 2 = X 2 {N). In turn, this allows us to hnd a simple solution to 
Eq. (2.11) in terms of N 


2 ^ ^ (2c|,(e/ - 3) - 3e/ -k 6) \ 

Af y 2c^{cTimp -k m2) j 


C 2 + 


F2{N) 

foFi 


(4.16) 


where C 2 is an integration constant, determined by the condition e = 1 at V* = 60, and 


Fi = [2c|, (e/ - 3) - 3e/ -k 6]^ , 


(4.17) 
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F2 (N) — 2c|) (c|) — l) \2(^rnp€f + 2fh? [A^ (ej — 3) ej + 3] — 3 (ej — 2) {Nej — 1)| , 

(4.18) 

Fg = m2 [24 (e/ - 3) - 3e/ + 6] . (4.19) 

In the following we find the range of parameters {cx>, m, ej) using the CMB constraints 
on {ns,r). Firstly, since the sound speed is constant we obtain the same constraint as in 
Sec. 4.1.3, namely c® — 0.98 to keep Ug within its observed range. 



Figure 7. Contour plots in the plane (m , e/). In the top panel, light and dark blue regions represent 
the 68% and 95% CL for the spectral index Ug, respectively. Black lines represent contours for different 
values of the tensor-to-scalar ratio, as indicated. In the bottom panel, the blue region depicts the 
95% CL for the spectral index Ug. We use c-p = 0.980. 

In Fig. 7 we depict the parameter space (m, e/) consistent with observations of the spec¬ 
tral index and tensor-to-scalar ratio. Taking cp = 0.98 and enforcing Ug = 0.968 ± 0.006, our 
plot shows that it is indeed feasible to obtain a tensor-to-scalar ratio as low as r ~ 6 x 10“^. 
Nevertheless, the plot also evidences that this requires a considerable tuning between m and 
ef. We have checked that using the 2a interval for the spectral index does not contribute to 
enlarge significantly the space where r ~ 10“^. In the absence of the aforementioned tuning, 
expected values correspond to the range 10“^ ^ ^ ^ 3 x 10“^. Moreover, we have checked as 
well that the space where r ~ 10“^ becomes incompatible with the observed spectral index 
even for small deviations away from cjy = 0.98. Consequently, finding r ~ 10“^ requires 
the combined tuning of m, ej and cp. Nevertheless, it seems fair to say that, despite these 
tunings, the DBIG model of inflation represents an improvement, albeit a moderate one, with 
respect to the DBI model studied in Sec. 4.1.1. 
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In addition, we verify the equilateral non-Gaussianity by using the approximate expres¬ 
sion for in Eq. (4.12). Since we consider a tiny variation of the warp factor we can 
practically neglect its contribution to non-Gaussianity. From Fig. 8 we can conclude that the 
DBIG model with varying warp factor leads to non-Gaussianities within the current obser¬ 
vational bounds. Gonsequently, we conclude that after including a varying warp factor the 
DBIG model of inflation could be of crucial importance with respect to B-mode detection 
and non-Gaussianities in future GMB experiments [7]. 



Figure 8. In this plot, we depict the non-Gaussian parameter as a function of rh (with m in 
units of mp). We take c-p = 0.98 and e/ ^ 10“"^ (Blue line) and e/ ^ 10“® (Green line). In this plot 
0.326 <a< 0.33 for 1 < m < 20. 

We finish this section by depicting the predictions of DBIG inflation for different sets of 
values of the model parameters in Fig. 9 and by summarizing our results in Table 1. We recall 
that the values collected in the table were obtained taking by enforcing the scalar spectral 
index to lie within its observed range = 0.968 ± 0.006 at the 95% GL and taking W = 60. 


Inflation 

r 

r/nt 

m^/mp 

GeV 

//rap 

DBI limit 

(0.01, 0.1) 

(-4.8,-0.7) 

6.63 X 10"® 

(0.95,1.82) 

~ 10^^ - 10^^ 

Galileon limit 

(0.13, 0.15) 

(-8.1,-7.93) 

< 0 

(0.64,0.70) 

~ 10® 

DBIG 

(0.01, 0.1) 

(-7.95,-7.5) 

< 0 

(1.7, 2.1) 

~ 10® - 10® 

(ii)Varying f 

(0.0068, 0.0095) 

(-7.95, -7.85) 

(2.41, 2.9) X 10"^ 

(5.9, 6.4) 

(6, 9) X 10^° 


(0.0018, 0.0027) 

(-8.01, -7.95) 

(3.6, 5.2) X 10~® 

(4.1, 4.6) 

(2.2, 3.4) X 10® 

(0.0006, 0.0007) 

(-7.63, -7.52) 

(1.52, 1.58) X 10“® 

(2.8, 2.85) 

(0.17, 0.18) X 10'^ 


Table 1. Inflationary observables in various limits of DBIG inflation. 


5 On a class of background solutions 

Until now, we have explored solutions to the background Eqs. (2.9) and (2.10) in which 
the sound speed and warp factor are either constants or time-dependent functions with very 
slow variation, although not simultaneously time-dependent. This choice is motivated by the 


- 15 






















Figure 9. Predictions of the DBIG model for iV* = 60 along with the Planck TT+low P+BKP+BAO 
constraints on the space {ns,r) at the 68% and 95% CL. The black line represents the case with 
constant sound speed and warp factor (c^) = 0.985, 1 < fh/mp < 1.25). Different model predictions 
for a constant sound speed and varying warp factor are plotted in red (c-p = 0.985, m = 15mp and 
5.1 < 10^6/ < 8.5), blue (cp = 0.98, fh = 15 top and 1.5 < lO^e/ < 2.6) and green (cp = 0.98, 
rh = 13mp and 0.07 < lO^e/ < 0.11). 


simplicity of the perturbation spectrum imprinted in the CMB, which strongly favors the 
simplest inflationary models. Nevertheless, it is reasonable to conjecture, and to some extent 
expected, that in the early stages of inflation, when the observable cosmological scales are still 
deep within the horizon, the background dynamics has been much different from the simple 
slow-roll evolution supported by CMB observations. Therefore, it is interesting to investigate 
what kind of inflationary dynamics does the DBI-Galileon model give rise to when the sound 
speed and warp factor become time-dependent functions simultaneously. In general, however, 
it is not possible to integrate the equations of motion for general functions cx>{t) and f{t). 
Owing to this difficulty, in order to find analytical solutions of the background equations 
we pursue a phenomenological approach in which we consider two different ansatze for the 
functions Ai and A 2 . 

If we allow the sound speed cp and warp factor / to change (ep, ej 7 ^ 0) the coefficients 
Ai ,2 become time-dependent functions. In such case, Eq. ( 2 . 11 ) can be rewritten as 


din H 
Ai - A2i/-2 


din a. 


(5.1) 


In what follows, we discuss two different parameterizations for Ai ,2 to find approximate solu¬ 
tions for a{t). 


Parametrization 1 

The simplest strategy to integrate Eq. (5.1) is to rewrite Ai ^2 as functions of H. Thus, we 
consider the temporal dependence for Aj (with i = 1 , 2 ) of the form 

A* = XiH ^^, (5.2) 

where Xi,ai are constants. Using this ansatz, Eq. (5.1) can be integrated to give 

2 F 1 ( 1,1 + p-,2 + I3\ ^ = A 2 (^2 - 2) In |Ka| , /? =- — -(5.3) 

V A 2 / a2 — ai — 2 
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where 2-^1 is the hypergeometric function and k is an arbitrary constant. Note that in the 
limit Q!i ^2 —)• 0 we can use the identity 2 .^i(l, l',2]z)z = — In |1 — 2 ;| to arrive at Eq. (2.15). 
Given the complexity of the above solution, substituting H = d/a to integrate the resulting 
differential equation in terms of a{t) is of no practical use. Thus, it is necessary to resort to 
numerical methods to integrate it. Nevertheless, if |/3| < 1 an approximation to the evolution 
equation is given by (see Appendix A for details) 


In 


1 - 


XiH^ 


X 2 


~ In 


Ka\ 


with 


^ (2 - a2)Ai 
“ 1 + /3 


(2 — a2)Xi 


(5.4) 


For 02 ^ Cl(l), the condition |/3| <C 1 implies |ai| <C 1. Provided H does not change 
exponentially, which can be certainly applied to the regular solution plotted in Fig. 1, we 
can approximate Ai by a constant since Ai ~ Ai (1 + ai ln(Pf/Pf*) + ...). This reasoning can 
be applied to the singular solution as well whenever it hnds itself sufficiently away from the 
singularity at t = t. Using Eq. (2.11), we rewrite Eq. (5.4) as 


j^2+ai a2 _ gign(Ai) (l + sign(ij)|Ka|"^') , ( 5 . 5 ) 

which can be integrated to obtain the scale factor a{t) in terms of hypergeometric functions. 
The implicit function (for simplicity we present the solution for k = 1 and vanishing ai) 
which dehnes the scale factor is given by. 


Ai (t — t) « — sign (h\ + sign (pf)) ( —sign (Ai) 


A 2 fa(2-“2)Ai + sign (H 


Ai 


2F1 ( 1,1; 1 + ;-sign 2 )Ai j ^ 02/2 


(5.6) 


From Eq. (5.5) we easily recover the background solution with constant sound speed and 
constant warp factor, Eq. (2.15), in the limit ai ^2 —^ 0. An important aspect of Eq. (5.5) 
is that it only requires |ai| to be small, whereas \a 2 \ can be relatively large, thus allowing 
a signihcant evolution of A 2 during inflation. Note that if we consider cj) constant, for 
consistency with the smallness of ai, then from Eq. (2.13) it follows that the evolution of 
A 2 is to be attributed to the warp factor /. Below we study the behaviour of the computed 
solution for different values of 02 . In view of Eq. (5.5), we may consider three cases consistent 
with > 0: 


Ai > 0 and H > 0. This case is illustrated in the left panel of Fig. 10, where for 0:2 < 2 
we have a singular solution when t —)• t . Any other solution with a 2 > 2 is regular at 
t = t. 


• Ai > 0 and H < 0. This regime takes place provided (|/t|a)^ < 1. A thorough numerical 
study of this scenario shows that only for a limited range of values of 02 the integration 
of Eq. (5.5) yields a well behaved physical solution for the scale factor. In the central 
panel of Fig. 10 we depict the solution for a few values of 02 in the range 3.5 < 02 < 5 

• Ai < 0 and H < 0. The constraint now is (|K|a)"^ > 1. This case, depicted in the right 
panel of Fig. 10, possesses smooth solutions for 02 > 2. Moreover, for large values of 
02 , the scale factor follows approximately a power law a (t) 
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Figure 10. Evolution of the scale factor a{t), according to Eq. (5.5), for Xi,H > 0 (left panel), for 
Ai > 0, < 0 (central panel) and for Ai,ii < 0 (right panel). For simplicity we take k = 1. 


Parametrization 2 

A second, simple alternative to solve Eq. (2.11) with time-dependent Ai ,2 is to parametrize 
their dependence as 

Ai = A,*(a/a*)"A (5.7) 

where A,*, are constants and a* is the scale factor when the largest cosmological scales exit 
the horizon. Defining 

z = lnH , y = ln(a/a*), (5-8) 

we find the exact solution to Eq. (2.11) (see Appendix A) 




2 A 2 

cti 


Q!l 


exp 


2Ai 

CHl 


,yai 


r I 

Oil q;i 


-|- Kexp 


— (e^"i - 1) 
ai 


(5.9) 


where r(s, x) is the incomplete Gamma function [46]. In the limit a \^2 —t 0 we easily recover 
Eq. (2.15), whereas for 0:1,2 7 ^ 0 we can use the asymptotic formula r(s,x) ~ when 

X S> 1. In such case, the above equation becomes 


g2z ^ I ai) _|_ ^gxp 

Ai* 


2A 


1* 


Ol 


(e^“i - 1) 


(5.10) 


If we focus on the background evolution while cosmological scales are exiting the horizon then 
0 < ?/ < 9, and yaj <C 1 provided |oi| <C 1. Neglecting higher orders in yoi we obtain 









(5.11) 


which can be integrated to obtain a{t) in terms of hypergeometric functions, and also gives 
Eq. (2.15) in the limit 01,2 —t 0. The implicit function (again, and for simplicity, we present 
the solution for k = 1 and vanishing ai) that determines the scale factor a{t) is given by 




2Ai^ (X 


— 0^2 / 

Ai. 


+ a 


2Ai 


2-ci 


(i' 


2Ai*— 0 : 2 ’ 4Ai*—2 q;2 


— 0:2 


+ 1 ;-- 


^2Al, 


A 2 * 


(5.12) 


a2A2* 


Similarly to Eq. (5.5), \a 2 \ is allowed to take on relatively large values in Eq. (5.11). 
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Notice that in Eq. (5.12) the hypergeometric function 2 F 1 is undefined when -\a 2 ^ 

is a negative integer. Therefore, from a formal point of view, by taking 02 < 0 we avoid the 
regions where Eq. (5.12) is undefined. In addition we must also impose that 02 7 ^ 2Ai*. When 
a 2 < 0, we have from Eq. (5.11) that H becomes singular if the scale factor a goes to zero. 
It can be checked in Eq. (5.12) that {t — t) is zero whenever a is zero, which amounts to have 
an undesirable singular solution for H when t = i. 

In view of our results, it seems reasonable to conclude that the solutions obtained using 
the ansatz in Eq. (5.2) for H < 0 (with either sign of Ai) provide a more appropriate qualita¬ 
tive evolution for a{t) than those described by the ansatz in Eq. (5.7). Therefore, our analysis 
demonstrates that, within the context of DBI-Galileon inflation, it is possible to envisage an 
early inflationary stage during which the warp factor undergoes a significant variation. The 
relevance of this result is that such phase can be smoothly connected to the last phase of 
slow-roll while allowing a marginal variation of the warp factor and agreeing with current 
CMB observations. In this sense, it is very suggestive to imagine that the early phase of 
rapidly evolving geometric structure could be connected to the very beginning of inflation. 

6 Conclusions 

In this paper we study the DBI-Galileon inflationary scenario, which constitutes a generic 
extension of the DBI model involving an induced gravity, and obtain the gravitational field 
equations allowing the sound speed cd and warp factor / that the model depends on to be 
time-dependent. We find exact solutions to the background Eqs. (2.9) and (2.10) when c® 
and / are constant. We obtain a singular behaviour at finite time for the scale factor and 
Hubble parameter when A 2 < and also a regular behaviour when A 2 > (see 

Fig. 1). We focus on inflationary scenarios under the slow-roll approximation and constrain 
the model parameters using the Planck 2015 results. In addition, we constrain the warp 
factor in the different inflationary regimes using GMB data. Notice that the warp factor scale 
might be important, regarding warped string phenomenology, to understand extra dimensions 
and warped geometries arising from string theory. We find that, in general, different warped 
geometries give rise to distinct inflationary predictions. In the case of constant cd and / (see 
Fig. 4), the tensor-to-scalar ratio is r > O (l0“^). Later, we considered the DBI-Galileon 
model with a slowly varying warp factor and find that the tensor-to-scalar ratio can be as low 
as r ~ 6 X 10“^ (see Figs. 7 and 9). However, we find that this requires the combined tuning 
of m, €f and cj). In any case, a varying warped geometry brings the predictions of the DBIG 
model closer to those of the Starobinsky model. 

Another aspect of our study is the violation of the standard consistency relation of 
single-field inflation, r = —Sn^. Since DBIG inflation is a class of generalized G-inflation, 
we find deviations away from the standard consistency relation r = —Sn*. However, with 
the exception of the DBI limit (see Fig. 2), the deviations found in the rest of cases under 
study are quite small (see Table 1). This result is consistent with the status about the tensor 
consistency relation in Galileon models as it is described in Ref. [40]. We emphasize that a 
prominent detection of the B modes, thanks to future GMB probes devised with a greater 
sensitivity, can discriminate models of DBIG inflation. 

Finally, we aim at describing an early stage of inflation taking place well before cosmolog¬ 
ical scales exit the horizon, we obtain general background solutions allowing an arbitrary time 
dependence for cx> and / by promoting the coefficients Ai and A 2 in the background Eq. (2.11) 
to time-dependent functions. To integrate the background equations analytically we pursue 
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a phenomenological approach, making nse of the ansatze in Eqs. (5.2) and (5.7). The validity 
of onr approximations demands that Ai remains approximately constant (ai ~ 0) for both 
ansatze, whereas A 2 can have snbstantial variation since 02 is not constrained to be small (see 
Fig. 10). This variation of A 2 , in tnrn, can be attribnted to a variation of the geometric warp 
factor / since cd remains approximately constant. From onr nnmerical exploration of the 
approximate solntion we conclnde that the ansatz in Fq. (5.2) provides a more appropriate, 
qnalitative evolntion for the scale factor. Onr analysis thns provides the intrigning possibility 
to consider an early stage of DBI-Galileon inflation (may be even connected to its very be¬ 
ginning) with a signihcantly varying geometric strnctnre that gives way, once the geometric 
strnctnre becomes approximately stabilized, to a hnal phase of slow-roll in perfect agreement 
with cnrrent CMB observations. 
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A Analytical approximations 


Parametrization 1 

Using the dehnition of the hypergeometric function [46] we have 

r (2 + /3) ^ r(l + n)r(l + P + n)z'^ 


2 F 1 (1) 1 + /3j 2 -|- /3; z) — 


E 


+ r (2 + f( + 6 ) „! 


(l + /?)x]f 


n=0 


+ ^ + n 


For /3 < 1 we can approximate 


1 


1 


1 - 


/3 


n + 1 — (5 

1 + 13 + n + 1-|-n 1-|-n ^ (n-|-1)^ ’ 

and snbstitnting in Eq. (A.l) we arrive at 


(A.l) 

(A.2) 


2 F 1 (1,1 + /3; 2 + /3; z) c. (1 + /3) f; (In |1 - z| + /? U^iz)) , (A.3) 

^ (n -I- ly z 

71=0 


where Li„( 2 :) = k~'^z^ is the polylogarithm fnnction [46]. Despite its being an excellent 

approximation for /3 < 1, snbstitnting the above into Eq. (5.3) leads to a differential eqnation 
still too complicated (to solve for a{t)) dne to the polylogarithmic fnnction L^iz). Onr aim, 
therefore, is to hnd a simple analytical solntion reprodncing the qnalitative behavionr of the 
scale factor. The simplest manner to achieve this is to neglect the term in the polylogarithm 
fnnction in Eq. (A.3). This simplihcation can be jnstihed after approximating 


00 
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71=0 


1 + y + n 
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In 11 — z 
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(A.4) 
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in Eq. (A.l), which holds provided /3 <C 1. In that case, after snbstitnting 2 ; —)> the 

resnlting backgronnd eqnation (Eq. (5.4)) has the advantage of being relatively simple. 

Parametrization 2 

Using the variables z and y defined in Eq. (5.8), onr Eq. (2.11) becomes 

^ q _ 5 ^ 

After mnltiplying by y{y) = exp [—(2Ai*/ai)e^“i], Eq. (A.5) becomes an exact differential 
eqnation 

df = P{y,z)dy + Q{y,z)dz = 0, (A.6) 

where 

P(y,z) =/i(y) - A 2 *e^“ 2 j ^nd Q{y, z) =-y{y) (A.7) 

Integral cnrves are of the form: /(y, z) = k, where k is a constant. Integrating / with respect 
to y in the first place we have 

/(y> = y 

where g{z) is to be compnted by demanding dzf{y, z) = Q{y, z). After integrating and solving 
for g{z) we find that the integral cnrves f{y, z) = k are determined by Eq. (5.9). 


B Perturbation spectra 
B.l Scalar spectrum 

To qnantify the amplitnde and tilt of the spectrnm we introdnce the variables dys = ^ dt, 
Zs = \/2a{J-sQsY^^ thanks to which the action Eq. (3.1) can be canonically 

normalized 

^ r r 7 " 1 

1 / . . , .0 .0 .. O 

00 for the 

(B.2) 


Sf"’ = n dys d^3 


(n')" - {VuY + 


Imposing the flat spacetime vacnnm solntion in the snbhorizon limit k/aH 
pertnrbation mode tt, we find we find 

u= "^V^sH^si-kys) , Us-^^Vsy- 

Using now = u/zg, we obtain the amplitnde in Eq. (3.5). 

To compnte the spectral index of the scalar pertnrbations 

Hg 1 = 3 2zzg, 

first we need to find i^g. Using the definition of we find 


(B.3) 




Ha 

Cg 


1 + 


fs + Qs 


+ ( l-e-^ + ^ 
' 2 2 


1 + 


fs + ds 


fsfP - OsQ?'’ 


(B.4) 
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The next step is to integrate dyg = {cs/a)dt. Assuming small and constant rj and rjg to 
neglect second order terms and integrating by parts we obtain 


Vs = 


1 + 


erj + esr]s 


{l-e-es)aH\ (e+ 65-1)2/’ 


(B.5) 


which is valid for somewhat large values of e and provided r], rjs are sufficiently small. Using 
now Eqs. (B.2)-(B.5) and neglecting second order terms in r],r]s one arrives at 


z 


1=(^\ 
Zs \ Cs J 


1 + 


fs + 9s 


1 — _ I 9s\ . f ^ . fs 9s 




( 2 ) 

S - 9s9s 


( 2 ) 


from which one can compute the spectral index in the linear approximation as 

46^ + 3 / 5 ^ ^ 5 ^ 

re.o - 1 - ■ 


(B. 6 ) 


(B.7) 


2 + 26^ + Qs* 

Using (3.6) we can compute the running index n( = Since we assumed ri,r]s 

approximately constant and small, the use of (B.5) allows us to write 


f VsaH dus 

n, = -- 


1 + 


ey + esr]s 


Cs din a (1 — 6 — 65 ) \ (e + 65 — 1 ) 2 / din a (1 — e — Cg) din a 


dus 


dus 


(B. 8 ) 


Provided r/ and r]s are small and approximately constant, we can expand (3.6) to first order 
in r],r]s. Using also Eqs. (3.4) and (B. 8 ) the running index becomes 




(<2\ 

(4 fs* + 9 s*) 2 ^ 9s* (1 + e*) 


(B.9) 


(2 — 2 e* — fs* + gs*)"^ 

B.2 Tensor spectrum 

Similarly to the case of the scalar spectrum, we introduce the variables dyt = ^ dt, zt = 
and Uij = Zthij so that the action in Eq. (3.7) can be canonically normalized 


•sf ^ = ^ / rfyt d^x 


K? - ( Vn ,,)2 + 


Zt 




Imposing the flat spacetime vacuum solution as in Eq. (B.2) we find 


Uij = Hjjl\-kyt) n 


2 _ 2^t 


(B.IO) 


(B.ll) 


where Cij is the polarization tensor. Using that hij = Uij jZt and taking into account the two 
polarization states, we arrive at the expression in Eq. (3.11). 


faH\ 


Zt I C 5 / 


ft + 9t 


+ 1 


-| + |-e + l ) + 


ft + 9t 


+ 1 - 


{u 


(2) „ „(2) 

t - 9t9t 


(B.12) 
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The spectral index of the tensor spectrnm is 

m = 3 - 2i't, 


(B.13) 


If slow-roll parameters are small, the hrst order approximation gives 

46^ T 3ft^ Qt* 

—2 + 2e* + /t* — gt^ 


(B.14) 
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